Study on X(3872) from effective field theory with pion exchange interaction 
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We study DD* {D* D) scattering in the framework of unitarized heavy meson chiral perturbation 
theory with pion exchange and a contact interaction. '^Si—'^Di mixing effects are taken into account. 
A loosely bound state X(3872), with the pole position being Mpoie = (3871.70 — i0.39)MeV, is found. 
The result is not sensitive to the strength of the contact interaction. Our calculation provides a 
theoretical confirmation of the existence of the 1"'"*' state X(3872). The light quark mass dependence 
of the pole position indicates it has a predominately DD* [D* D) molecular nature. When the n 
mass is larger than 142 MeV, the pole disappears which makes impossible the lattice simulation of 
this state at large quark mass. 



The narrow resonance structure named X(3872), dis- 
covered by the Belle Collaboration in the B+ ^■ 
ii'"'" J/^TT+TT^ process ill and then confirmed by CDF 
and DO Collaborations through its inclusive produc- 
tion in proton-antiproton collisions [2, l3|, has inspired 
heated discussions both experimentally and theoretically. 
In 2006, the Belle collaboration studied the B^ — > 
D^ D^ TT^ K^ decay process and found an enhancement of 
the D°D°Tr° signal just above the D°D*° threshold li], 
with the resonance is peaked at 



+0.3 . 



Mx = 3875.2 ± 0.71^;^ ± O.SMeV 



(1) 



A later analysis of Belle data including new data on 
D* ->■ D-f^ gave 



Mx = 3872. etoi ± 0.4MeV 



(2) 



The latest PDG value for the X(3872) mass from the 
J/'^X decay mode is [6| 



Mx = 3871.68 ±0.17MeV 



(3) 



The angular distributions and correlations of the 
Tr^TT~J/ip final state studied by the CDF collaboration 
indicate two possible quantum numbers of this state, 
J-^*-^ = 1++ or 2 '". The radiative decay reported by 
Belle and BaBar bolsters the 1++ assignment, while the 
37r invariant mass distribution in J/ipuj decays slightly 
favors 2 '". 

The line shapes of B+ ^ XK+ in the J/^tt+t:^ 
and D^ D^ TT^ / D^ D*^ modes and the corresponding pole 
structures have been studied independently by two 
groups [71, |8| . The effect of energy resolution effect was 
taken into account in [9| and a twin-pole structure was 
found, suggesting that the X(3872) can be identified as 
a 2^Pi cc state strongly distorted by coupled channel 
effects. The effects on the lineshape of X(3872) from 
nonzero decay width of D*^ and inelastic channels was 
taken into account by Braaten [10|. 

The proximity of the X(3872) to the threshold of 
D°D*° strongly suggests that the X(3872) isprobably 
a loosely bound D^D*° molecular state |lll - [l3 |. Other 



interpretations include normal charmonium [15l4l7| . a 
tetraquark state [l8| and a ccg hybrid |li| . In the above 
investigations, many kinds of potential models were used, 
leading to different conclusions. Recently, one-pion ex- 
change as a possible binding mechanism in the X(3872) 
was revisited 20] . The authors argued that it is not suf- 
ficiently binding for this purpose, suggesting other short- 
range interactions should be included and may be respon- 
sible for the X(3872) formation. 

Most of the theoretical studies can only provide 
a model and parameter dependent prediction for the 
X(3872). Until now the theoretical situation concerning 
the existence of X(3872) and its J^'-^ number was un- 
clear. Even for lattice simulations, there are also many 
lattice studies predicting 1++ X(3872) mass ranging from 
3850 to 4060 McV [21 24], but with various uncertainties 
of their own, where the key difficulty is the challenging 
task of extracting the excited states. The lattice sim- 
ulation indicates a lower mass with J^*^ number 2 ^ 

2AM- 



Effective field theory is a very powerful tool to study 
hadron properties at low energy. A pionless effective 
field theory describing DD* scattering was proposed 
in I2.6| , where the s-channel bubble diagrams are summed 
by solving Lippmann-Schwinger equations to produce a 
bound state. However, to obtain a pole near thresh- 
old, the low energy constants must be huge. In other 
words, for a small "natural" interaction, there is no 1++ 
X(3872) appears. The pion exchange interaction is well 
known from the D* to Dtt decay. We will see that with 
pion interaction alone, there is a bound state pole near 
threshold which is not sensitive to the strength of the 
contact interaction. 

In this Letter, we provide a reliable, parameter free 
calculation using unitarized chiral perturbation theory. 
This approach has been successfully applied in the light 
and strange quark sector [27| . It has also been applied to 
the heavy meson case with open charm [28|, |29[ . The key 
advantage of this method is that the Lagrangian is well 
defined. The mass and width of the resonance or bound 
state are obtained from the pole analysis of the scattering 
amplitude. One need not assume the constituent compo- 



nents of the resonance as in potential models. 

Here, we calculate DD* {D*D) scattering amplitudes 
up to one loop in the framework of heavy meson chiral 
perturbation theory. We will focus on the J^'^ = 1++ 
channel. The effects of D* finite width and ^Si-^Di 
mixing can be naturally taken into account. Pade ap- 
proximation is used to construct the physically unitarized 
T matrix. 

We first construct C-parity even initial and final states, 
which are a superposition of D*D and DD* , 



\X+>=^{\D*D>+\DD*>) 



(4) 



We shall look for possible bound state, virtual state or 
resonant poles of the following transition amplitude. 



r_ 



++ ^ {X+\T\X+) 



1 



{Tii+Tr2 + T2i+T22) , (5) 



where Tij are the four relevant transition amplitudes. 
The general Lagrangian describing 4-boson contact in- 
teractions can be written as [261 



£(0) 
C2 



p(Q)tp(Q)v-(Q)t^(Q)M _^ p(0)tp(Q)v-(Q)t^(Q)M 



p(Q)tp(Q)^(Q)t^(Q)M ^_ p(Q)tp(Q)y(Q)t^(Q)M 



(6) 



where P^^^ = {D",D+,D+) and V^'^'> 
(L)*°, £)*+,£)*+) are the heavy meson fields, while 
piQ) = {D°,D-,D-) and V'^^^ = {D*° ,D*- ,D*-) are 
the heavy antimeson fields. 

In order to include the effect from pion exchange 
on the properties of X(3872), one needs to introduce 
the covariant chiral lagrangian describing interactions 
between heavy mesons and Goldstone bosons, which can 
be written as (sil-Isl 



£« = 2.9.(T/W)tp(«)+p(Q)tyW)X 



2.9.(K^,^'tp(0)+p(Q)t^(Q) 



"ab^ 



where 



Uf^ = i{u^d^u - udfj_u^) 



exp( 



^0_ 
V2f' 



(7) 



(8) 



with (f) is the 3x3 matrix containing the Goldstone 
boson fields. F is the Goldstone boson decay constant 
in the chiral limit, which we identify with the pion 
decay constant, F = 92.4MeV. The coupling constant 
f/TT depends on heavy meson masses, g^r — g^/MuMij* , 
where the dimensionless constant g can be determined 
from the strong decay D*^ — > D^n^. Taking the rate 
from PDG @, we find g = 0.30 ± 0.03. 







(a) 



(b) 



FIG. 1: Tree level Feynman diagrams. Double lines indicate 
the vector D*" or I)*" mesons, solid lines pseudoscalar D" or 
I)" mesons. The dashed line represents the pion. 



The tree level scattering amplitudes from contact in- 
teractions have already been given in Ref. 26|. Tree 
level amplitudes from one-pion exchange should be added 
to the contact contributions, since both of them are of 
0(/), 



Tn 



-Cie*(p3)-e(pi) 



Ti2 = C2e*{p3) ■ e{pi) + 
T21 = C2e*{p3)-e{p,) + 



^gl £*(P3) ■P2e{pi) -Pi 
F2 u- M2 

4g^ e*(P3) ■P2e{pi) -Pi 



p2 



Ml 



T22 
where 



-Cie*(p3)-e(pi). 



(9) 



u{s,cos9) = {pi~piY 



Af2 



Md'Md 



^ 2p2cos0 



(10) 

with A = Md* —Md and p the external three momentum 
in the center-of-mass frame. 

In the following, we only focus on partial waves with 
J — 1 and omit this index for simplicity. The projections 
of tree level amplitudes at 0{p^) with I = I' = and with 
positive C-parity are given by 



T 



(0) 
-+.SS 



-(C2-C1) 



M 

3F2 



^^l 



2p 



,2 H^ 






(11) 

where /ij = rn^ — A^ is the effective pion mass. The 
superscript denotes the chiral order. In this work, we 
count p, Utt ^ C(p)- 

One loop diagrams with contact and D* Dn vertices are 
shown in Fig. [2l all of which are 0{p^) in a naive power 
counting scheme. In the following, we use a similar pro- 
cedure to that used in the analysis of NN scattering [35| 
to deal with the one loop diagrams. We will find that 
diagrams (a), (b), (c) and (d) start to contribute already 
at O (p) . The partial wave amplitudes have the usual uni- 
tarity cut along the positive real axis in the complex p^ 
plane, starting at p^ — 0, which is related to the s-channel 
D*D intermediate states. We finally get 



T 



(l),a 
+ + ,SS 



iiCi - C2f 



P 



87t{Ml 



M, 



D* 



(12) 
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=v 



(d) 
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(f) 



=r 






Al 




(9) 



(h) 



FIG. 2: One loop Feynman diagrams. The notations are the 
same as Fig. [T] 



The contribution to the elastic scattering anrphtude from 
the triangle one loop diagram with three relativistic prop- 
agators is shown in Fig. [^^b). Finally, this yields partial 
wave amplitudes with positive C-parity: 



T 



i),b 



++,ss 



= r' 



(l),c 



++,SS 



iC2 - Cr)gl 



[2tp-filToip)] 

(13) 

The contribution to the elastic scattering amplitude from 
planar box diagram is shown as Fig. [2][d). The other 
diagrams in Fig. 2, starting to contribute at Oip^), are 
neglected in our numerical calculation. 

The partial wave amplitudes satisfy the perturbative 
unitarity condition 



ImT 



(1) 



T 



(0) P 



'Stt. 



-.T 



(0)* 



++ 



(14) 



where T^l and T^^ are the '^S'l— ■^£'1 mixed amplitudes 
from which wc can construct the physical, unitary am- 
plitude using Fade approximation. 



T 



phy 



T 



(0) 



\T 



(0) 



T^] 



■T 



(0) 



(15) 



With our conventions, the relationship between S matrix 
and T matrix is given by 



1 



P 



nphy 



A-Ky/S 



Trnp) 



(16) 



The physical masses of the scattering particles 
are taken from PDG [6], Md = 1864.91, Md* = 
2006.98, TO^o = 134.98(in MeV). 

We first set X — C2 ~ Ci =0, finding that pion ex- 
change alone is strong enough to form a bound state 
just below D*D threshold. The pole position in the 
complex p-plane and corresponding pole mass are p = 



(-15.46 -hi24.62)MeV,M = (3871.70 - i0.39)MeV, re- 
spectively. Regardless of the small imaginary part, which 
is due to the finite width of its constituent, T{D* — > Dtt), 
our result is in good agreement with Ref. [ll[. That the 
obtained mass is very close to the experimental mass of 
the X(3872) is remarkable since we did not adjust any 
parameter. The only parameter g is fixed by the decay 
width of D* — > Dtt. On comparison in Ref. [26[, in order 
to get the correct mass of the X(3872), the particular 
value of \r (8.4 x 10""' MeV~^) is chosen, corresponding 
to our dimensionless parameter A = XrMuMo* = 3144. 
The obtained mass is very sensitive to A and the bind- 
ing energy is proportional to 1/A^. The binding energy 
decreases with the increasing interaction strength A. Es- 
pecially at A = 0, the binding energy is infinite. This 
behavior is not physically reasonable. Most important, 
due to the particular choice of A, it is not convincing 
that the X(3872) should exist. In our work, we want to 
determine the existence of X(3872) theoretically without 
using any information concerning the X(3872). With the 
pion exchange interaction included, the pole is automati- 
cally generated and the mass is close to the experimental 
data without adjusting any parameter. 

We now check the effect of including additional at- 
traction of shorter range by increasing the low energy 
constant A. The pole trajectory is shown in Fig. [31 from 
which we can see that the binding energy changes slightly 
as we increase the strength of the attractive contact inter- 
action. For example, the mass changes from 3871.70MeV 
to 3871.55 MeV when A increases from to 3000. It is 
clear that the A dependence of the pole is highly sup- 
pressed once pion exchange is included. It also shows 
that the pion exchange interaction is the main reason for 
the system to be bound. Without the pion exchange in- 
teraction, there is no pole around the DD* threshold with 
only a contact interaction except for the particular value 
of A. Our calculation theoretically confirms the weakly 
bound state X(3872) with J^^ number 1++ exists. This 
is a prediction directly from the unitarized heavy meson 
chiral perturbation theory rather than fitting the exper- 
imental data. 

Compared with potential models, the advantage of this 
pole analysis approach is that it does not specify the com- 
ponents of the state. The obtained state is a physical 
state which could be a mixing of many pure states with 
the same J^ numbers. In other words, we have no in- 
formation of the structure of the state. For example, for 
X(3872), we do not know whether the obtained state is 
a traditional meson or tetraquark. Therefore, we need 
other methods to obtain the information concernin g it s 



37|, 



structure. As emphasized in a series of papers [36 
the quark mass dependence can provide important infor- 
mation on the structure of a state. If the X(3872) is a 
pure cc state with no constituent light quarks. Its light 
quark mass dependence only comes from sea quark con- 



tributions, which should be very weak, as for the case of 
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FIG. 3: Pole trajectory as increasing the strength of contact 
interaction. 



Z?s (1968) shown in lattice simulations [38| . 

We fix the c quark mass at its physical value and vary 
the light quark masses (pion mass). At tree level, the 
pion mass dependence of D^*' mesons can be expressed 
as H^ 



MdM (M^) = MjjM \phy 



2hQ + hi 

MdM \phy 



{Ml-M%ny) 



where the values of /iq and hi are taken from [29|. The 
pion mass dependence of the mass of X(3872) as well as 
the threshold is shown in Fig. 2] From the figure, one can 
see that for a small change of quark mass (from zero to 
several MeV), there is an obvious change of the X(3872) 
mass which shows the light quark component in the state. 
In general, the mass of a loosely bound A—B molecular 
state is given by 



Ml, = Ma 



Mb - Eh 



(18) 



where Eb is the binding energy. The pion mass depen- 
dence of Eb is expected to be much weaker than that of 
Ma and Mb- Thus, the mass of an A — i? molecular 
state is almost the same as that of the A — B threshold. 
From Fig. 21 one can see that the pion mass dependence 
of the X(3872) is in good agreement with the expecta- 
tion for a. D* — D molecular state. We can also see that 
the tiny width decreases with increasing pion mass, since 
the phase space for D* decay to Dtt will be suppressed. 
Eventually, the bound state pole will be absorbed by the 
D*D threshold as /x^ — )• 0. 

When ttItt is large than A (142 MeV), there is no bound 
state or resonance pole. This is an important result re- 
lated to the lattice simulation. Recently, a lattice simula- 
tion result showed that there ig no l^"*" state near the ex- 
perimental mass of X(3872) 24*]. The authors explained 
that if X(3872) was a molecule state, it was not expected 
to be found in their simulation. Our result shows if the 
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FIG. 4: Upper plot: Pion mass dependence of the mass of 
X(3872) (solid line) and D* D threshold (dashed line); Lower 
plot: Pion mass dependence of rx(3872)/2. 



simulation is at large quark mass, it certainly can not be 
found for any l"*""*" operator, even though the operator 
has a large coupling with the molecule state. 

In summary, we studied the DD* {D*D) scattering 
using unitarized heavy meson chiral perturbation theory. 
The standard one pion exchange is included in the chiral 
Lagrangian. The parameter g is determined from the 
PDG value for the decay D* -^ Dtt. The obtained pole 
mass is very close to the experimental data which is not 
sensitive to the low energy constant A. This confirmation 
of the existence of the l'''^ X(3872) is a parameter-free 
prediction rather than a fitting of the experimental data. 
The light quark (pion) mass dependent behavior shows it 
is probably a D* — D molecule. The disappearance of the 
pole at large pion mass makes it impossible for current 
lattice simulations to find the 1++ X(3872) . Future 
lattice simulations must be carried out at the physical 
pion mass in order to find it. 
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